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EXPANSION PROBLEMS OF ORDINARY LINEAR DIFFERENTIAL 

EQUATIONS WITH AUXILIARY CONDITIONS AT 

MORE THAN TWO POINTS* 

BY 

CHARLES E. WILDER 

1. Introduction 

The boundary value and expansion problems for the equation of the nth 
order with boundary conditions at two points have been studied by Birkhoff. f 
BocherJ has suggested the generalization of these results to the equation 
with auxiliary conditions at more than two points. Such generalization of 
the essential properties of the differential system has been carried out by the 
author, and in this paper is given the proof of the convergence of the expansion, 
which may be studied quite independently of the other results. The formal 
development of the boundary problem and a more detailed discussion of the 
form of the series will be presented in other papers. 

The differential equation is taken in the form 

d n u d n ~ 2 u 

-^ + * + P2 "^ZJ + • • • + fn U + \U = 0, 

in which fi, f%, • • • , p n , are functions of the real variable x , continuous 
together with their derivatives of all orders in the closed interval ( a , b ) . 
The auxiliary conditions apply to k points of the interval (a, b ) , of which 
the first and the last are the ends of the interval. These points, arranged in 
order of ascending algebraic magnitude, will be denoted by 

a\ = a , a 2 , a 3 , • • • , a^-i ,a k = b. 



* Presented to the Society, April 29, 1916. 

f Birkhoff, Boundary value and expansion problems of ordinary linear differential equations, 
these Transactions, vol. 9 (1908), pp. 373-395. See also Tamarkine, Sur quelques 
points de la theorie des Equations diff&rentielles liniaires ordinaires et sur la generalisation de la 
serie de Fourier, Rendiconti del Circolo Matematico di Palermo, 
vol. 34 (1912), pp. 345-382; and Birkhoff, Note on the expansion problems of ordinary linear 
differential equations, Rendiconti del Circolo Matematico di Palermo, 
vol. 36 (1913), pp. 115-126. 

% B6cher, Boundary problems in one dimension, Proceedings of the Fifth International Con- 
gress of Mathematicians, Cambridge, 1913,vol. 1, pp. 163-195; see p. 166. 

Trans. Am. Math. Soc. 87 415 
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The auxiliary conditions may then be written 

Wi{u) = W u (u) + W u {u) + ••• +W ki (u) =0 (i = 1,2, •••,«), 
in which 

Wh (u) = a 3i w (ay) + a'ji «' (a y ) + a^ u" (ay) + • • • + a^r 1 ' « ( " _1) (ay) , 

and ayi, %,-, ayi, • • ■ , ctfr ]) , are constants, real or complex. 

It is possible to define explicitly in terms of the auxiliary conditions and a 
fundamental system of solutions of the differential equation, a function 
6 (x, s; X) , having the properties of a Green's function of the system. These 
properties will be discussed elsewhere with the exception of a single one, 
which forms the subject of this paper, namely that 

(1) f(x)=\ha±j'JG(x,*-,\)f(*)d8ik, 

where T is a contour in the X-plane enclosing the first r poles of G (x, s; X) , 
and / ( x ) satisfies certain conditions of continuity, in some cases of the sort 
usually imposed, in others more restrictive, according to the nature of the 
auxiliary conditions. The auxiliary conditions themselves are subject to 
restrictions beyond those found necessary in the two point case. This integral 
may be regarded as the sum of the first terms of an infinite series of residues, 
corresponding to the series of orthogonal functions in the familiar cases, and 
these residues can in many cases be expressed linearly in terms of the char- 
acteristic functions of the differential system. Since the integral will be con- 
sidered here on its own merits this interpretation will not be considered 
further. 

2. Preliminary formulas 

A fundamental system of solutions of the differential equation may be 
chosen analytic in X. Let ?/i(x,X), ?/ 2 (x,X), •••, i/ re (x,X), be such a 
system. Then let D (X) be the determinant which has Wi{y,) for the ele- 
ment in the zth row and jih column. Its roots are the characteristic numbers 
for the system, but will concern us only as points where an expression having 
D (X) in its denominator ceases to be analytic. Let Di (s , X) be the deter- 
minant which has the (n — i )th derivative of t/y ( s ) for the element in the 
ith row and jth column, and let Z>2 (x, s; X) be this determinant with the 
elements in the first row replaced by the yj(x) with corresponding subscript. 
Then let 

(2) G{x,a;\) = ±lD t /Di, +ifz>5, -Ux<s. 

If now iV(x,s;X) is the determinant formed by bounding 2)(X) on the 
top by the yj(x) with corresponding subscript, and on the right by the W% ( G ) 
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with corresponding subscript, and filling in the upper right-hand corner with 
G(x, s; A), then the function G(x, s;X) mentioned above is given by the 
formula 

(3) G(x,s;\) = (- l) n N/D. 

The proof of the convergence of the integral is based on the use of the 
following asymptotic solutions of the differential equation, due to Birkhoff.* 
Make the substitution X = p" . Then on any region S of the p-plane, defined 
by the relation 

It (l+l)ir 

- S arg p S , 

in which I is a positive integer less than n > 1, there exist n linearly inde- 
pendent solutions, analytic in p , and such that 

Vi{x,p) =u i (x,p)+e"^- a ^l, 

dyi dut „,_ „, Ei 

(4) dx dx ^ p™' 1 ' (i-1,2, •••,!»), 



dx"' 1 (fa*- 1 ~ l ~ p m-n+l> 

where 



Ui(x,p) = e" w ' (x ~ a '> ■ 



1 T + ' ' " + 



p p" 

Here m is any positive integer, the w's are the roots of the equation 

w n + 1 = 0, 

the functions «,,- ( x ) are finite and continuous for all values of x in the interval 
(a, b) , and £ is a generic notation for functions of p (and other variables) 
uniformly bounded for large values of p . In part of the work it will be neces- 
sary to use a more explicit form of the functions «•,-,■ (x) , due to Milne, f viz., 
Uij(x) = <j>j(x)/(wi) j . The functions <j> are then the same for all the n 
solutions ?/i ( x , X ) . 

By the substitution X = p" the X-plane goes over into two adjacent S 
regions, S% and S 2 . The integral (1) becomes 

2r~if f n P n ~ 1G ( x ' s >P")f( s ) dsd P 



+ 2 



■"■■ f f np^G(x,s;p n )f(s)dsdp, 



* Birkhoff, loc. cit., p. 381. 

f Milne, Note on asymptotic expressions in the theory of linear differential equations, Bulle- 
tin of the American Mathematical Society, vol. 23 (1917), pp. 166-169. 
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in which 71 and y 2 may be taken as arcs of the same circle with center at the 
origin. 

Before substituting these solutions in the auxiliary conditions, the condi- 
tions are normalized as follows: 

First, the number of conditions of order n — 1 at the points a and b is re- 
duced, by means of linear combinations, to two at most; then, the remaining 
conditions of order n — 2 at these points are reduced to two at most; and 
so on. 

It is convenient to introduce Birkhoff's notation, by which [w] represents 
an expression of the form 

a b c E 

P P P 3 P 

in which w , a, b , etc., are independent of p but need not be constants, and E 
is a function of p and other variables bounded for large values of p . 

The result of substituting the asymptotic solutions into the auxiliary con- 
ditions is 

Wu(yj) = (pwy )*"[««], 

W 2i (yj) = (pw 3 y*[a 2i ]e^- a >, 

********** J 

W ki (yj) = {pw j ) h «\a ki ]e<' a >«'- a \ 

Because of the normalization of the auxiliary conditions hu = hu = ha, 
j = 2 , 3 , • • • , k — 1 . So let us set hu = hu = h , and we have A,_i is hi 
S hi+i , and no three successive h/s are equal. 

In the odd order case let n = 2p, — 1 . Then the indices 1 to n can be 
determined in such a way that for any p in S 

R(pWi) ^ R(pw 2 ) Si R(pwz) ^ • • • =1 R ( piv n ) , 

where R(w) indicates the real part of w . Hence for S 



R(piv. : ) < 


(» = 1,2, ••■ > /i-l), 


>0 


(1=0 + 1,0+2, ••-,«), 


= 


(i = and arg p = -~ 2 *- J , 




(i-1,2, •••,0-1), 




(.7=0 + 1,0 + 2, •■-,«), 



and it follows that 

Wi(yj) = {piOjf'Wu] 

= {pw i )H<x k i]e> w > <b - u) 

= ( pw s ) h ' [ a u ] + ( pwj )"' [ a ki ] e"**- + £ ( P™i ) h " I an 1 e^"'""' 

0=m). 
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In the even order ease let n = 2p. Then the indices 1 ton can be deter- 
mined in such a way that for any p in S 

R(pwi) < (1 = 1,2,3, ■■-,/* — 1), 

>0 (i= M + 2, ■■■,n), 

= (i = n,ii + l, and either arg p = — or arg p = x ) , 

and since — w M = w, + i , we have 

Wi{ yj ) = (pwy )**[««] (i = i,2, ...,„- i), 

= ( pw,- )* [ a ki ] tr**-* (j = m + 2 , • • • , n ) , 

= ( p«v ) fe [ a u ] + ( pw M )*■ [ a ki ] e""M<*-«> + Z ( P«V )*« [ a K ] e<"V«"- a > 

1=2 

O'-m), 

= (-/», )'" [««] + (- P«V ) fe [ «H ] e -"%(^°) 

£-1 

+ E ( - p«v )*" [ «k ] rf - ""* (ai_a) ( y = *• + 1 ) • 

3. The function D (X) 

Preparatory to taking up the main convergence problem it is necessary to 
study the distribution of those points where the function D ( X ) vanishes or 
becomes nearly equal to zero. This study is simplified by considering separ- 
ately the cases where n is odd and where n is even. 

Odd order case. When n is odd, the substitution of the asymptotic ex- 
pressions for the W's in the determinant D (X) gives 

D(X) = f{p h < f[ e'^^A, 

in which A is that determinant which has for the element in the iih row and 
jth column, [ «h ] wf for j = 1,2, • • • , p — 1 ; [ a*,- ] wf f or j = p. + 1 , p. + 2 , 
• • • , n , and finally 

B.- = «£([a u ] + [« K ]A (M + £ (pwj^'la,;]^^'), 

i=2 

for j = p . Let y be the highest power to which pw^ enters in any Bi ; then 
A can be expanded in the form 

(6) A = [B 1 ] + [d k ] e<"V^°> + ( pwv )- J g [ 0, ] e"'V (««-«) j . 

We shall now prove the 

Theorem. 7/ zn (6) neither di nor 6k is zero, and if for v =% 1 neither 62 
nor 0jt_i is zero, then D (X) Aas an infinite number of roots, and there exist two 
absolute constants w and £2 such that the number of roots of D (X) in any circle 
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about the origin of radius r is less than fir and greater than ur and, finally, for 
values of X uniformly away from the roots, D (X ) remains uniformly away from 
zero. 

This last statement is made more explicit by inequalities (16) and (21) below. 

We assume, then, that neither 0i nor 0k is zero. The number v may be 
negative, zero, or positive. If it is negative A can be written more simply 

[0i] + [0 k ]e pw ^ a K 

Birkhoff* has shown that the roots of this are asymptotically equal to 



'-(-I) 



2iic4^~i 

w^ib — a) ~^ w u (b — a) ' 

and that A remains uniformly away from zero if p is uniformly away from a 
root. 

When v is zero, A becomes 

(7) [Oi] + [B k ] <?<"»* <*>-«> + £ [ 6i } e ,m ^ a ^ a) , 

1=2 

and now at least one other in addition to 0\ and 0k is assumed to be different 
from zero. In treating this case the distribution of roots of the function 
obtained by suppressing the brackets is determined, and then the roots of A 
are proved asymptotically equal to these. 

It is convenient to make a change of variable. Let pw M ( a 2 — a ) = z , 
and divide through by 0i . The equation to be considered is, then, 

(8) f(z) m 1 + Q e* + Qi e°* + Q 2 e°* + • • • + Qk-2 e°*-*> = 0. 

The Q's are complex constants, and in particular Qk-2 and at least one other 
are different from zero. The c's are real constants such that 1 < Ci < c 2 
< • • • < cjb-2 . Let z = x + yi, and let R ( <j> ) denote the real part of <j> , 
and I ( <p ) the coefficient of the pure imaginary part. Let e be an arbitrary 
small positive quantity, and choose X so large that forx = — X ,f(z) differs 
from 1 by less than e, and at the same time for x = X , f{z)/e Ck ' % " differs 
from Qk-2 by less than e. The region between x = X and x = — X will 
be called the region R throughout this paper. We shall suppose that e is 
chosen small enough so that f(z) does not vanish outside of R . 

Now f(z) is analytic in any finite region. Hence in any interval, no 
matter how small, along the y-axis, it is possible to find a line y = constant 
along which /(z) does not vanish. Let y = Y x and y = Y 2 be any two such 
lines. The number of roots of/ (z) in the rectangle x = X ,y = Y\,x = —X, 
y = Y 2 is determined by following the change in arc/(z) as z describes this 

* Birkhoff, loc. cit., pp. 383-5. 
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rectangle. Now 

arc/ (2) = sin \fr> 

and for y constant I (/) has the form 

I (/) » d e* + dt e CiX + d 2 e c **+ ■■■ + 4-2 e°™* , 

in which the d's are real constants. The finite roots of I (/) are the same 
as those of 

h (/) = do + 4 e ( *-» x + 4 e«*-» x + ■ ■ • + 4-2 e^- 1)x . 
But between any two roots of this is a root of the derivative 

IJ (/) ■* 4 (ci - 1 ) e^" 1 ^ + 4 (c 2 - 1 ) e^* 

+ ••• +d^ 2 (c fc _2-l)e (c *- a - 1) % 
while its roots are the same as those of 

h (/) - 4 ( ci - 1 ) + & ( c 2 - 1 ) «<*-»" + • • • + 4-2 ( c^ 2 - 1 ) e««-** . 

The argument is repeated till finally Ih-z is reached, which has only one term 
and does not vanish for a; finite. Thus I (/) has at most (k — 2 ) finite roots. 
Should any of, the Z's have multiple roots the argument must be slightly 
changed, but the result is, of course, the same. 

Now as z traces the line x = — X , /(z) remains very nearly equal to 1 , 
so that the resulting change in arc / (z ) is less than some small positive quantity 
i?i which approaches zero with e . As 2 traces the line x = X the change in 
arc/(z) differs but slightly from c*_2 ( Yi — Y\) , as is at once apparent 
from the formula 

f(z) 
arc/(z) = arc e°*-* + arc"^. 

In other words, the change in arc along this line can be written 

Ch-i{Yi - Yx) ±V2, 

where 172 is a small positive quantity which approaches zero with e. Finally 
along the lines y = Y\ and y = F 2 the function sin {arc/(z)} vanishes at 
most k — 2 times, whence arc/(z) increases by less than (k — l')ir + t\z 
and (k — 1 ) ir + rjt respectively, where 173 and j?4 are again small positive 
quantities which approach zero with e. To sum up, as z traces the rectangle 
the resulting change in arc / ( z ) is not greater than 

c^2(F 2 - Y 1 )+2(k-l)ir + r l 
nor less than 

Ck-i(Y 2 — Yi) — 2(k — l)7r — i\ (i = >n + 7)2 + 7)8 + ';4), 

if Yi — Y\ is large enough to make this expression positive. Since 77 can be 
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made arbitrarily small by taking e small enough, it can be dropped from these 
inequalities. So finally in any interval of R of length I the number of roots of 
f(z) is not greater than c k -2 l/2-ir + k — 1, nor less than c k - 2 l/2ir — k + 1 , 
if this number is positive. In particular, by making I very small, one sees 
that a multiple root of / ( z ) cannot be of order higher than k — 1 . 

It is now necessary to prove that if z remains uniformly away from a root 
of / ( z ) then / ( z ) remains uniformly away from zero. This is proved by 
means of the following auxiliary theorem : 

Theorem. Given f(z, xi, x 2 , • • • , x n ) , continuous in all of its arguments 
and analytic in z for a t Si x ; Si b{ and z in a closed finite region S of the complex 
plane, and such that for no set of values of the x's is the number of roots of f(z) 
greater than a given constant N; if for any set of values of the x's, z is a point at a 
distance greater tlian 5 from a zero of f(z) and from the boundary of S , then 

\f(z,x 1 ,x 2 , ■■■ ,x n )\^\, 

where X is a real positive constant independent of the x's. 

It is no restriction to assume that 5 is small enough so that for- each set of 
values of the x's the corresponding region of variation for 2 , defined by the 
hypotheses, actually contains points. This follows from the assumption 
that for any set of values of the x's the number of roots of / ( 2 ) is less than iV . 
Now for any particular set of values of the .t's there is a X , such that for these 
x's and for 2 restricted to its corresponding region of variation, 

|/(2, xi, xi, ••• ,x n )|S X, 

while there is a value of 2 for which the equality holds. We have thus defined 
a function of the x's, X ( Xi , x 2 , •••,»„), which, by means of the easily proved 
fact that the zeros of f{z, Xi, x 2 , • • • , x n ) are continuous functions of the 
x's, can be shown to be continuous in the closed region a; si Xi si &,• . Since 
X ( Xi , x 2 , • • • , x n ) is continuous in a closed region, it has a lower limit which 
it attains. Call this limit X. Then for some value of z and the x's, \f\ = X. 
Hence X is not equal to zero, for / does not vanish in any of the regions for 
which the function of which X is the lower limit was defined. Hence the 
theorem is proved. 

First apply this theorem to the function 

f(z, <t>u 02, • • • , 0*-2) = 1 + Qo e° + Qi e^ + ^ + ■■■ +Q k -2 e c ^ + *-*. 

The region S can be that portion of R for which Si y Si 2ir , augmented by a 
strip of width 8 surrounding it. Obviously the breadth of R, that is, the 
value of X , depends only on the Q's and not on the <f>'s. It is also clear that 
the function satisfies the conditions of the theorem. Let X indicate the 
number determined by the theorem corresponding to S . 
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Consider, then, the function (8) in R. Any interval of R, 

2fo^y^2(l + l)r, 

can, by a transformation y — y + 2lir be brought to the interval ^ y ^ 2w , 
and / becomes 

f(z + 2ln) = 1 + Q e z + Qi e^' 11 ^ + ■■■ + Q k - 2 J 1 -*-* 1 " -* 1 . 

But, without changing /, each expression 2lircji can be replaced by a <j>ji 
where 0S<A = 27r. So / ( z + 2lwi ) has the form above, and for z at a 
greater distance from a root than S , \f ( z + 2lir ) | == X . By giving I all 
integral values, all of R is covered by these intervals, and so finally, for z in R 
and at a distance from a root of / ( z ) greater than 5, |/(z)|=£X>0. 

We can now prove that the roots of (7) are asymptotically equal to those 
of (8). First it is clear that by choosing e small enough we can be sure that 
all the roots of (7) lie within the region R . Let us denote (7) by \{/ ( p ) and 
(8)by/(p). Then in R 

-A(p) ^ e l f( P )+- p , 

from which it is obvious that the roots of ^ ( p ) are very nearly equal to those 
of / ( p ) for large values of p . Now 



(9) 
whence 



arciHp) = arc 0i + arc/(p) + arc f 1 + f~fT~r~ ) • 



Let S be an arbitrary small positive quantity. Then according to the 
previously determined distribution of the roots of / ( p ) , it is possible to 
enclose these roots in groups of less than k in simply connected regions of 
major diameter less than fcS , and such that no point of the boundary is at a 
distance from a root of / ( p ) less than 5 . On the boundaries of these regions 
we have then |/ ( p ) | = X , where X is determined by the theorem above. 
Hence for large values of p , the expression E/( di f ( p ) p ) is very small on these 
boundaries, and so the resulting change in arc (1 + E/( di f ( p ) p ) ) as z traces 
the boundary of a region is zero, and the change in arc \p ( p ) is the same as 
the change in arc / ( p ) . So the number of roots of \p ( p ) in any such region 
is for large values of p the same as the number of roots of /(p) in it; that is, 
the roots oi\j/(p) are asymptotically equal to those of / ( p ) . 

From this together with (9) can be deduced the fact that for large values 
of p in R the function \j/ ( p ) remains uniformly away from zero when p is 
uniformly away from a root of \p ( p ) . Outside of R the function ip ( p ) of 
course remains uniformly away from zero. 
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It is convenient to consider next the case in which A reduces to 

[0i\ + [0k] e" w ^ b - a) + (pwj'[e a ] e > V«*-> , 

which, in particular, it always does if k = 3 . It is assumed that none of the 
0's are zero, and that v =£ 1 . Then 

A = [02KV«*- a >Ao, 
where 

(io) A ^ ( P w M y + [ e b ] ^v 6 -* + [e a ] e -pv*-«> , 

and in this 

It is necessary to consider separately the part of 8 in which R ( pw M ) is positive 
and the part in which it is negative. 

R ( pw M )g0. It will now be shown that in that part of S where R ( pic,, ) 
Si , the roots of Ao are asymptotically equal to those of 

/(p) = ( P « v :r + a <r<"v o *-»>. 

To find the roots of /(p) set pw^ = x + yi = re'*, and — O = r a e'*". 
When / ( p ) is equated to zero the second term may be transposed, and then, 
when absolute values and angles are equated, the two following equations 
result: 

r* = r a e*< a ^ > , v <fi = $ a - y(a 2 - a) + 2mir (m = 0, + 1, ±2, ± 3, •••)• 

The zeros of / ( p ) are the points of intersection of these two curves. Since, 
in that part of the p plane under consideration, x S 0, y > 0, the first curve 
is equivalent to 

y = VrI"e- 2 *< a *-°>/" -x\ 

and this has but one branch which is in character much like the exponential 

curve, and which for large values of p approaches parallelism to the line 

R(pw„) = 0, i. e., 

lim (dx/dy) = 0. 

p=00 

It is also easily shown that 

lim (x/yY = 0, 

which is to say that for points on this curve 

lim <f> = x/2 +2lir (I =0, + 1, ±2, ± 3, •■•)• 

p=00 

When this is substituted in the second curve it is found that the roots of 
f ( p ) are asymptotically equal to points whose y coordinates are given by 

_ <j> a + 2mir — v (|tt + 2lw ) 
y ~ o 2 — a ' 
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and since v is an integer this is equivalent to 

<f> a + 2mw — hvw 

V = . 

at — a 

So the roots of/(p) lie along the curve 

(11) |(pwJl = |e a e-"V— )| 

at intervals whose lengths approach 2ir/(a 2 — a) as p becomes infinite. For 
large values of p these roots are all simple, for it is obvious that / ( p ) and /'( p ) 
cannot both vanish for the same value of p when p is large. 

In order to prove that in that part of /S where R ( pw>„ ) = the roots of A 
are asymptotically equal to those of / ( p ) , it is necessary to show that |/ ( p ) | 
remains uniformly away from zero when p is uniformly away from a root. 
For this purpose let T be the region included between the curves defined by 
the relation 

|(pw / ,)1 = Ke~'" v(a2 ~ 0)±,r l- 
Let U be the region between T and the line i?(pw M ) = 0, and let V be the 
remainder of S for which .R (pw„) < 0. Choose a as any convenient positive 
real number. The following inequalities are immediate: 

for U , | ( pw„ )' | sS e" \ 6 a tH"**^ \ , 

for V , | ( pw„ y | S e'" 1 6 a e-cV^-) | ; 

from which it follows that 

fortf, |/(p)|S(e*-l)|0 ae ->V— >|, or (1 - e-')\{ P w,y\, 
(12) 

forF, |/(p)|^(l - g-*)\0 a e-"**-*\. 

For p in T it is necessary to proceed differently. Let p be any root of 
/(p) =0, and call those points of T which are at least as near to po as to 
any other root of / ( p ) = , the domain of p . Now the strip T is asymp- 
totically of constant width, and it is easily shown that this width is 2a /( a%—a). 
On the other hand the distance between two roots of /(p) =0 was shown to 
be asymptotically equal to 2ir/(a 2 — a), so that finally if po is far enough 
out I he distance of any point of the domain of po from po will be less than 
2 {a + 7r)/(a 2 — a) . In this domain/(p) can be written 

/ ( p ) ^ / ( p ) -/ ( po ) - ( pw, y - ( po «v )' + 0» e-'V«*-«> - e a *-«v*^> 
(13) = < (p - po) (p- 1 + p"- 2 po + • • • + pr 1 ) 

+ d a e-'V '*-'*) ( 1 - gO-nOVi— ) ) . 
In T the absolute value of the first term is less than 



2(o- + 7r) 
j, 

02 — a 



, 2(<T + 7T) 
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and for sufficiently large values of p this is less than a constant times | po | r ~ 1 . 
Moreover, since in T , | pw„ | " is between two constant multiples of | e~ l ""^ a2r ~ a) \ , 
it follows that the first term is less than a constant times 

I g-[(i'-l)/>']p"' (/ .(.«2-a) I # 

In the parenthesis in the last term let 

(p — p )w M = r( cos </> + i sin 4> ) , where — ir Si</> Si w . 
The parenthesis becomes 

/ 1 „r(az-a)(eos<l>+ism $) \ 

which can vanish only at the points given by 

2/tt 

cos (6 = 0, sin <f> = — -, ; 

r (a 2 — a) 

and the point p . In the domain of p there is only one of these points, namely 
Po itself, and for that portion of the domain whose points are at a distance 
not less than 8 > from p the absolute value of the parenthesis as a function 
of r and <f> is continuous and positive in a closed region. Hence it has a 
positive minimum which may be denoted by r\ . The last term of (13) is then 
greater than a constant times 

I g-P«V( a 2- a ) [ 

for that portion of the domain of p at a distance from p greater than S. 
This, combined with the inequality on the first term, gives the result that 
for large values of p in T and at a distance from a root of / ( p ) = greater 
than or equal to 5 > , 
(14) \f(p)\>\e-">^-^\K 0> 

where K o is a positive constant. Similar inequalities have been proved for 
any p in U or V , so that this, with a new constant K i , may be assumed to 
hold for all that part of S for which R ( pw^ ) Si f or large values of p at a 
distance from a root greater than or equal to 5 . 
From (10) is obtained the relation 

Ao = / ( p ) + - e -» V*- a > + [ b ] e PV 6 -<"3> , 
P 
which can be factored into 



/ E fT' V' 2 - a > [ 9 b ] e" X' | - a2 > \ 

\ 1+ 7 SiP) + '~Kp) ~)' 



(15) A =/(p)^l + 

Consider the portions of T , U , and V , bounded by two circular arcs \p\ = p' 
and | p | = p" , (p' < p") , such that no point of either arc is at a distance 
less than 8 from a root of / ( p ) = . Call such a region S . It will be shown, 
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by considering the change in arc A as p traces the boundary of S , that the 
number of roots of Ao and of / ( p ) in S is the same, provided the innermost 
boundary of S is taken far enough out. Consider first the last term of the 
parenthesis in (15). Throughout the region S the numerator is bounded. 
In U the denominator, f(p), can be replaced by (1 — e~°) (piv^)" from 
inequality (12), and so for large values of p in U this last term is small. From 
the second of the inequalities (12) it can similarly be shown that this last 
term is also small for large values of p in V , while from inequality (14) the 
same result is obtained for the boundary lines of S where they cross T , pro- 
vided the innermost boundary of S is far enough out. On the other hand, 
from inequalities (12) and (14) the expression e~ pw ^" r '' I) /f(p) is bounded on 
the boundary of S , and when it is multiplied by E/p the result is small when 
the innermost boundary of S is far enough out. Under these conditions, 
then, the parenthesis in (15) can be written as 1 + 77, where 77 is small enough 
so that the arc of 1 + V will return to its original value when p traces the 
boundary of S . It follows that the change in arc A is the same as the change 
in arc / ( p ) , so that A and / ( p ) have the same number of roots in S . 

Now exactly the same result is obtained if p traces the circumferences of 
circles of arbitrarily small radius about the roots of / ( p ) , and this, combined 
with the previous result, shows that in that part of the region S for which 
jR ( p«v ) = > the roots of A are asymptotically equal to those of / ( p ) . 

An important inequality for A can be deduced from formula (15). For p 
large and at a distance from a root greater than 5 inequality (14) holds. For 
T and V the exponential can be replaced by a constant times (pw^Y , while 
for U the resulting inequality for/(p) has already been obtained in (12). 
The parenthesis in (15) can be made nearly equal to 1 , and in particular, 
greater than 1/2, by taking p large enough; so that the inequality is obtained: 

|A | >|pwv|" K 2 , 

in which K 2 is a constant depending only on 5 . The result of applying this 

to A is (see (10)) 

(16) |A|>|(ptO"^V«3-«>|Z 3 , 

in which K3 is a constant depending only on 5. This holds for p in that 
part of S where R ( pw M ) is less than or equal to zero and greater than some 
fixed p , and at a distance from a root of / ( p ) at least as great as 5 . Equally 
important is the result, easily deducible from the same formulas, that under 
the same conditions A is greater than a constant depending only on 5 . 

R ( pwv ) = . By entirely similar reasoning it is found that for that part 
of S where R ( pw^ ) IS the roots of A are asymptotically equal to those of 



(ptO' + fcC-M 



,pw„(6-a») 
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which lie along the curve 

at intervals whose lengths approach 2w/(b — a 2 ), and which are all simple 
for large values of p . Inequalities on / ( p ) , A , and A , similar to those above 
are found also in this case. 

Finally the odd order case with k points is in general, i. e., when certain 
constants do not vanish, merely a combination of the cases already con- 
sidered. The characteristic equation A = is as given in (6). The coef- 
ficient of (pw M )' can be brought to the form (8) by dividing by e'V 8- " 5 
and setting pw„ ( o 3 — a% ) = z . It is now assumed that none of the con- 
stants 6%, 02, 0*-i, Ok, is zero. It follows then, that the roots of this coef- 
ficient are all in a certain region R lying along the line R ( p«v ) = , and that 
for p at a distance from a root of this coefficient not less than S the coefficient 
is in absolute value greater than | e^" 2- ") | K± , in which Ki is a constant 
dependent only on 5. For convenience let this coefficient be denoted by Aj. 
Then the inequality is 

(17) |Ai|>|e p V a «-)|X 4 , 

and A is 

A = [0x] + [0*]e<"V«—> + (ptO'Ai. 



This can be written 

(18) A = (p« V yA 1 (l + 



[0i] + [0fc]e p V^ a) ' 
(pwJ"A x 



Now in R the exponentials are bounded. Hence, if in R any region be drawn, 
such that no part of the boundary is at a distance from a root of Ai less than 
5 > 0, then on the boundary Ai remains uniformly away from zero, and if 
the region be far enough out the parenthesis in (18) has, on the boundary, 
the form 1 + V, where \y\ < 1 . From this can be deduced the fact that the 
roots of A in R are asymptotically equal to those of Ai . 

Outside of R for that part of S in which R ( pw^ ) < the roots of A are 
asymptotically equal to those of 

/(p)^(p^r + |e- (, v«-") = o, 

and for that part of S in which R ( pw„ ) > they are asymptotically equal to 
those of 






(P«V )» + ■— e'V* 1 -"") = 0. 



To prove the first of these statements, write 

A , [0i] 0i , , [0*le' , V 6 -°) 
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[0] - 6iT,[dj]e<"°»( a >-°* 



(19) =/(p) + 



^=3 



[fl^e'V*-^ 



Ai 2 Ai 



*-i 



np; v + Aifl,/( P ) Ai«,/(p) + A t /( P ) ;• 

The roots of /(p) have already been found, for it is the/(p) of (11). To 
prove the roots of A asymptotically equal to those of / ( p ) , it is obviously 
sufficient to prove that the parenthesis above has the form 1 + i\ , for p on 
the boundary of any region outside of R and in that part of S where R ( pw^ ) 
< , and such that the distance of any point of the boundary from a root of 
/ ( p ) is not less than S > , provided that this region is taken far enough out 
from the origin. For such values of p the inequalities (12), (14), and (17) 
hold, and they will be used to show that the parenthesis has that form. From 
(14) and (17) is obtained |/(p)Ai| > K^Ki, and since the symbol [0] is 
the same thing as E/p it is evident that the term [0]/Ai dzf(p) is small for p 
large. In the next term the same inequality on / ( p ) Ai is used for T and V , 
since in these regions, for p large, E(pw M ) is large and negative, and the 
numerator, and so the whole term, is small. For U the inequality |/ ( p ) | 
i£ (1 — e~")\(pw li y\ from (12) is used for /(p). The remainder of the 
term may be written in the form 

fig— P» M {(°2-«)— (°S-<>!!)} _ 2£g— f>« M (os-«) [1— («8-"2)/(«2-«)]_ 

Now in U the inequality 
holds, from which it follows that 

I e _ pt V (aS_a) t 1 "" ( ° 3-< ' 2) /( a 2-")3 I < I ( pyj yll-( a 3~<*!!) /(02-«)]'| J£ 

where K is an absolute constant. The exponent of this power of p is cer- 
tainly less than v, so that because of the factor Kpw^)"] from /(p) in the 
denominator, the term approaches zero as p becomes infinite. By reason of 
inequality (17) the last term of the parenthesis can obviously be written 
E/f(p). In T and V inequality (16) shows that this is small for large values 
of p, while the inequality | /(p) | = (1 — e~")\{pw ll ,) v \ from (12) serves 
the same purpose for U. Thus the parenthesis has been shown to be of the 
required form for the values of p in question, and so the roots of A are 
asymptotically equal to those of / ( p ) for that part of S outside of R for which 
Ripw,,,) < 0. 
The proof for the remaining part of S on the other side of R is similar. 
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It is to be noted that the above reasoning applies to the whole of that part 
of S in which R ( pw^ )£0, even including the strip R , provided that p is 
at a distance greater than 5 from a root of /(p) or Ai, and is large enough. 
Hence from (19) can be obtained the following inequality which holds for all 
p's in absolute value greater than a certain p and at a distance from a root of A 
greater than 5: 

(20) \A\>\A 1 \-\f( P )\K 6 , 

in which Ka, is a constant depending only on S . 

In particular this shows, by reason of inequalities (14) and (17), that if in 
this region p remains uniformly away from a root, then A remains uniformly 
away from zero. 

Another inequality for A can easily be deduced from the above. First 
replace |Ai| by |e"V^- a )|Z 4 from inequality (17). Then for the region U 
replace |/ ( p ) | by ( 1 — e~" ) | ( pw M ) | " from (12), and for the regions T and V 
replace it by | e -' > V a ^- a ) j K\ from (14). But in the regions T and V the 
exponential can be replaced in this inequality by a constant times | ( pw„ ) | ", 
so that finally the following is obtained : 

(21) \A\>\(pw ll Ye^ a ^\K e , 

in which K$ is a constant depending only on S , and the inequality holds for 
all values of p greater in absolute value than a certain p and at a distance 
from a root of A greater than 8 . 

Even order case. When n is even the substitution of the asymptotic ex- 
pressions for the W's in the determinant D (X) gives 

Z)(X) ^IIV* ft e'^^A, 

4=1, j=n+2 

in which A is that determinant which has for the element in the ith row and 
the jth column, [ au ] wf for j «= 1,2, • • • , p. — 1 , and [ olm ] wf for j = p, + 2 , 
;u + 3 , ••• ,n, and 

Di^v^ilaxi] + [« K ]^^» + £ (pw ll ) h »- h >[a K ]e'"°»^) y-M), 
(21a) '= 2 w 

F i = (-wJ & ([« li ] + [« w ] e -"V^) + X:(-p^) fa '- fe [ce i; ]e- p V<"-a)) 

where n = 2p . 

It is evident from the definition of the bracket symbol that in general this 
determinant expression for A can be expanded in the form 

(22) A - Pt ( pw, Y* + P 2 ( pw.Y + P 3 , 
where the highest and lowest order terms of the P's are 
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Pi: [0ii]e"V a *-i- a 3> and [e li )e-" w ^ a ^- ai >; 

P 2 : [ 6>2i ] e"V^) and [ 8 22 ] e'""^"-^ , if b - a 2 > a k _i - a , 
(23) 

or [0 31 ]e<"V a »-''> and [ 3 2 ] ^"V""-") , if b-a 2 <a k ^-a; 

P 3 : [04i]e p "> ( ^> and [ 42 ] r"V Hl » . 

The number p is the exponent of the highest power of p appearing in a Dj 
or Hi . 

We shall now prove the 

Theorem. If none of the d's in (23) are zero, then D (X) has an infinite 
number of roots, and there exist two absolute constants co and Q such that the 
number of roots of D (X) in any circle about the origin of radius r is less than 
flr and greater than ur , and finally for values of X uniformly away from its roots, 
D (X) remains uniformly away from zero. 

This last statement is given more explicitly by inequality (25). 

Should the number v be negative, A reduces formally to Birkhoff's even 
order case,* and the problem of distribution of roots will not be considered 
farther here. If v is zero, A becomes merely a sum of exponentials which 
when multiplied through by the non-vanishing factor e 9W ^ b ~ a ' > all have real 
parts of the same sign, and so are reduced to the form (8) already treated. 
In either of these cases the formal expansion (22) is not needed. But when v 
is positive we make use of this expansion, which tacitly assumes that at least 
two of the Di's and H/s actually involve this power of p . This is, of course, 
equivalent to assuming that the highest order derivative in the auxiliary 
conditions at points other than the end points appears in at least two con- 
ditions. We also assume that none of the d's in (23) are zero, which is to 
say that certain of the determinants from the matrix of coefficients of the 
auxiliary conditions must not vanish. The method here used can be applied 
equally well with very slight modifications to many of the cases thus excluded; 
but a detailed study of the more specialized cases is beyond the scope of this 
paper. 

The more general even order case, then, is treated as follows: The function 
Pi can be reduced to the type (8) by multiplying through by e p V a *~ 1-02) and 
making the proper substitution. Its roots, then, all lie in a region R , and the 
proof that the roots of A in R are asymptotically equal to those of Pi is im- 
mediate. To find the roots outside of R write 

The line R ( pw^ ) = is now not the bisecting ray of each of the 8 regions 

* Birkhoff, loc. cit., p. 387. 
Trans. Am. Math. Soc. 28 
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as in the odd order case, but the line of demarcation between two of them. 
Hence the region R overlaps the t wo S regions into which the X-plane is carried 
by the transformation X = p" . Consider first the region Si in which R ( pw^ ) 
< . The above expression can be written in the form 



Tr{w + m'~"-* 



L "32 J 



+ (p«v )'& + &• 

Suppose that a 2 — a is not equal to b — ak-i and, for convenience, let a 2 — a 
> b — ajb_i . Then the term of highest order in Qi is 

[6] e~i' w ^ b ~ a '" l) , 
and in Q 2 , 

[0]e-"V a »- a >. 

The expressions Qi and Q 2 are merely those terms which must be added to 
the product of the two factors in order to make the result equal to A/Pi . 

The roots of the two parentheses are found by the methods already de- 
veloped, and it can be shown by means of inequalities of the type (12) and 
(14) that the roots of A in Si and outside of R are asymptotically equal to 
those of the two parentheses. 

Similarly the roots of A in S 2 and outside of R are asymptotically equal to 
those of two parentheses of the same type, but with positive instead of nega- 
tive exponents. 

The following inequality is easily found for large values of p which are in 
Si and at a distance from a root of A greater than 8 : 

|a| > i(p W ,r + [^ e -pv—)| . [(pw.y +|^-j^v»— >| . \Pi\k,, 

in which Ki is a constant depending only on 8 . If 

A = Ae pM " (&-o) , 
the inequality for A is 

|A| > I Pi e'V" w -"'» I • I (pw>„)" e'V* 2 -") 

(24) +[Si-i^'-' l —+ral K '- 

It can also be proved that for large values of p in Si which are at a distance 
greater than S from a root of A , the absolute value of A remains greater than 
a positive constant depending only on 8 . 

When o 2 — a = b — ak-i the two systems of roots lie along two curves of 
the type (11) in which both have the same exponential, although in general 
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the coefficients are different. If these coefficients should have the same 
absolute value the roots all lie along one curve, and if they are identical the 
roots are all double. An inequality like (24) holds in this case also, the only 
difference being in the 0's and the value of K 7 . From (24) it is easy to deduce 
the inequalities 

I A | > | ( pw„ ) 2 " e" V«*-°> e" V 6 -*") | K s , 

(25) >\{pw ll )»e» w * ( - a ^ a '>\K<>, 

>\(pw li Y e^^- tt ^\K w , 

which hold under the same conditions as (24), and in which K s , Kg, and K w 
depend only on 5 . 

4. The expansion problem 
The expansion theorem is as follows : 

Theorem. Given f(x), any real function which together with its first v 
derivatives (or if v < 1 with its first derivative) is continuous in the closed interval 
(a, b) , then its expansion in terms of the residues of the Green's function (3) 
converges to f(x) at every interior point of (a,b) , if the auxiliary conditions 
have the following properties : 
In the odd order case 

Neither 0i nor Ok in (6) is zero. 

If v i£ 1 then in the coefficient of (pw^)" neither 2 nor dk-i is zero. 
In the even order case 

None of the 6's in (23) are zero. 

The intervals b — a 2 and Uk-i — a are longer than any other of the subinter- 
vals of (a, b) included between any two of the points a , a 2 , a 3 , • • • , au-i , b . 
The essential purpose of this last condition is merely that the second point 
from either end shall be farther away from that end than the first point from 
the other end is from that end. The conditions on the 6's require merely that 
certain determinants from the matrix of the constants in the auxiliary con- 
ditions do not vanish. The condition on the 0's in the even order case is 
stated in the above form in order to have it as simple as possible. If v < 1 
it is not necessary that all these 0's be different from zero. 

To prove the theorem it is necessary to consider separately the odd and 
even order cases. But before that is done it is convenient to introduce a 
new form of the expression for G(x, s;X) given in (2). If G be expanded 
in terms of the first row of the numerator determinant it can be written 



in which 



— 1 " 

G(x,s;\) s= ±~J^yi(x)yi(s), 



(26) *(.)- npn - x [ 1 +^ W -+l^+---+Z pW ^ + p -n)- 
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This expression for the y's in so far as it differs from the one given by Birk- 
hoff * is due to Milne, f In the work to follow, if v < 1 , Birkhoff's form would 
be sufficient, and the convergence proof is no more complicated than in the 
two point case. The complications necessary in case v = 1, are due to the 
necessity of using the more explicit formulas of Milne. 
The above form of G gives for W% ( ) 



W t {Q)m -lY,Wu(yj)y } (s) ±l£,W ii {yi)y j (s) ± ■■■ 

£ j= -\ 4 j=i 

1 n 1 H 

±,E Wh- U {yj)yj (s)+;E Wm (y,)vj (*) , 

in which the plus sign is taken for values of s to the left of the point at which 
the term in question applies, and the negative sign for values of s to the right. 
Odd order case. R ( pw„ ) si . The convergence of the integral (5) is 
now to be studied. It has already been separated into the sum of two inte- 
grals, one in Si and one in S 2 . In the odd order case it is necessary to separate 
each still farther into the sum of two integrals, one in that part of S where 
R(pWp) < 0, and the other in that part where R(pw li ) > 0. The suc- 
cessive paths of integration are arcs of circles | p | = constant extending across 
both the regions Si and S 2 , and so chosen that no point of such an arc is at a 
distance less than a preassigned 5 > from a root of D . That such arcs 
can be drawn in such a way that the sum of the orders of the roots between 
two successive arcs is not greater than 2k when 5 is chosen small enough, 
follows from the distribution of the roots of A as determined earlier in this 
paper. The parts of these arcs in Si and S 2 will be denoted by 71 and 72 as 
in (5). Let S represent either of these regions. Then let B be the point 
where the line R ( pw„ ) = in S cuts one of these arcs, and let A be the other 
end of the arc that crosses the region where R ( pw„ ) < , and C the other 
end of the arc that crosses the region where R ( pw M ) > . The convergence 
of the integral over the arc AB will be studied first; that is, the integral 



(27) 



= 2nj f n P n ~ l °( x > s >P n )f( s ) dsd P- 



Before making use of the expression (3) for G ( x, s; p n ) it is convenient to 
alter the last column by adding to it each of the columns 1 to p multiplied by 
the term ^y, (s) with corresponding index, and each of the columns p. + 1 to 
n multiplied by the term — \yi ( s ) with corresponding index. The resulting 
last column is: 

1st element, 11yi(x)yi(s) for y < x, or — ]C y% (%)yi(s) for y > x; 



* Birkhoff, loe. cit., p. 391. 
t Milne, loc. cit., p. 167. 
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(1 + i)th element, 
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- E Wu(yj)yj(s) + 

(28) j =k +1 



- E W2i(yj)yj(s) 
J2W 2i (y 3 )yj(s) 



y + 



+ 



- E W k -i, i (yj)y ] (s) 
JlWk-i,i(yj)yj(s) 



+ 22W ki (y i )y ] (s) T 



j=i 



in which the upper term of each brace is taken for values of 5 to the right of 
the point to which that term applies, and the lower for values of 5 to the left. 
The asymptotic solutions are now substituted in the expression , and then 
the factor p hi of D is divided into the (1 + i)th row of the numerator, and 
the factor e pw>(6_a) into the jth column. The factor — 1 is divided into the 
last column. Thus is obtained 

in which N is the determinant A bounded on the top by e l>v "^ x ~ a) [ 1 ] with 
corresponding index for the columns 1 to /x, and by e pm(x ~ b) [ 1 ] with corre- 
sponding index for the columns p + 1 to n. On the right the ith row is 
bounded by Ft and the upper right-hand corner is filled in with F , where F, is 



(29) '=* +1 '= 2 



and F is 



- E [anvfr+^ePA*- 



j'=n+i 



+ T,aiiU'j" +1 ]e (,u ' ia '- s) 



j=i 



+ Da b -^ +1 ]e« M , 



j=i 



£eP»i(— > [ w .] > for s < x, and - E e' w ' ( *- s) [wy], for « > x; 

and for convenience we give Bi, the ith element in the ;uth column of A: B t is 

'"*' ( [ a u 1 + E ( P w » ) hH ~ hi I «k ] e" V--> + [ ttw ] e oV»-> ) . 



w„ 



This expression for G (x, s; p n ) can now be substituted into (27), and since 
s occurs only in the last column of N , the function / ( s ) can be multiplied 
into this column, and the s integration performed. It is then to be noted 
that all terms of the determinant save possibly those in the ;uth and last 
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columns are bounded. Moreover the cof actor of F is — A , and so the deter- 
minant can be expanded to give 

i ft <* =1 n n r r h \ 

+ ~. D EE^^ M, x *■,/(*)& Up (c) 



(30) 



It will be remembered that £ is a generic notation for functions of p (and 
other variables) uniformly bounded for large values of p. In particular, 
in (6) and (c), Ej r i is identically zero for r = 1 . 

The integral (a) will be treated first. It is necessary to separate it into the 
sum of two integrals, one from o to a; and the other from x to b . The result 
of integrating these by parts with respect to s gives (n/4n)f(x) + e and 
((p — l)/4n)/(a;) + « respectively, for, as is easily seen, 



/ 

JAB 



p a e" a ' ix - a) dp , for j<n, 

>AB 

(31) and 



/ 

J Ah 



p a er v *- V) dp , for i>p, 



approach zero as the radius of the arc AB becomes infinite, for a any real 
constant, and x distinct from a and b . One also easily shows that 



(32) ——dp 

Jab P 



approaches zero as this radius becomes infinite. 

To evaluate (b) and (c) it is only necessary to notice that 1/A is bounded, 
while both D r and the s integral are certainly less than a constant times p" . 
Hence (6) and (c) reduce merely to bounded functions multiplied into integrals 
of the type (31), and so approach zero as the radius of AB becomes infinite. 

The evaluation of (d) is more complicated. We now make use of the 
more explicit asymptotic solutions of Milne, and his expressions for the y's, (26). 
In the formula for F t in (29), the square bracket [a u wj" +1 ], for example, 
has the more explicit form 
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Consider now 

f Fif(*)da. 

Ja 

It can be broken up into the sum of several integrals, since Fi is a sum. Of 
these integrals consider first the one with the exponential involving ai ( I =}= 1 , 
k ) . When the more explicit form of the bracket symbol given above is 
used this becomes 

p H-M ( f'j^au wf' +1 e" w ^ a '- e) f (s)d* 
\ J a y=i 

+ I - E auwf i+1 e< M ' <ai -' ) f(s)ds (a) 

Ja y=l r=l (, pWj ) 

(33) 

+ f- t n tf~4^^w^f(s)ds (6) 

Ja, j=n+l r=l (. pWj ) 

+ rt,§;e<'W-*->w? +1 f( S )d S 

J a J=l P 

- E ~ n e< m ^ l - >) w h / +1 f(s)ds). (c) 

a, i=i*+l P / 

Now integrate (a) by parts and add the results, obtaining 

— ( / ( m ) E ( - wf ) + £ wj« / ( a ) e p^( a «- ffl ) + £ «$» / ( b ) ««<«*-») 

P \ i=i i=i i=n+i 

E w)«e"^ < "- 5 > f'(s)ds+ I - E «J" e<"">< ar -*> /'(«)<&). 

Since n > &k > 0, the coefficient of /(a;) is zero. The remaining integrals 
can again be integrated by parts,* and again the terms involving ai will drop 
out. This is repeated hu — hi times, which will give the factor au/p h "~ h ' 
outside, and this is just enough to cancel the factor in (33). Then absolute 
values are taken, /c A »- ; ''- 1 ) (s) is replaced by the upper limit to its absolute 
value, and the exponentials are again integrated, giving E/p. The total 
result may be written 



Eh 
P ' 

Next consider the rth term of (b). The factor l/p r can be removed from 



(34) E En p*"-**- 1 gp^-") + Y, En p''"-*'- 1 e^e*-*) + — . 

i-l j=H+l P 



* It is to be noted that the number* of derivatives here required has been provided for in 
the hypothesis of the expansion theorem, and also that if v < 1 all this work with Milne's 
formulas is unnecessary. 
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under the sign of integration, and au r (s)f(s) can be considered as one 
function. The integral can then be treated as was (a), but it requires only 
hu — hi — r integrations by parts and one further integration. The result 
can then be written in the same form as (34), though in this case, to be sure, 
the E's are not only bounded but actually approach zero. The terms (c) are 
obviously of the form E/p, so that the whole of (33) can be written in the 
form (34). Finally a single integration by parts reduces the terms involving 
the points a and b to the form E/p , so that the integral 

Fif(s)ds 
can be written as 



r 



(35) Z ( E Em p""-"^ 1 e^ '-) + E E in p*"-*'- 1 «**««-»> ) + — . 

1=2 \j=l ^'=H+1 / P 

When this is substituted in (30d) the last term here gives a result of the form 
(32). The other terms for j 4= ju give results of the type (31) with the value 
ai for x , and since for these terms ai is distinct from a and b , the result of the p 
integration approaches zero as the radius of AB becomes infinite. The only 
terms left from (35) then are 

*— i 

y ft phi-M-l gP«V(«!-a) 
1=2 

which are multiplied by {E/i^)e pv '^ x ~ a ' ) for the p integration. If it can be 
proved that 

*— 1 ft. „*tf-*t-l gPw^ai-a) 

S £ 

is of the form E/p, then this last term will give a result of the form (32), 
and (30d) will have been shown to approach zero as the radius of AB becomes 
infinite. Now inequalities (16) and (21) prove this for v > 0, while for 
iiSO the power of p is negative while the exponential is bounded and A re- 
mains uniformly away from zero, and the same result is obtained. 
To sum up the results, for the radius of AB large 

The work for the case R(pw a ) > runs along parallel to this in every 
detail, and the result is obtained that for the radius of BC large 

he = i/(z) + e- 

This result holds for both Si and S 2 , so that (1) converges to / ( x ) in the 
odd order case. 
Even order case. In the even order case the paths of integration are the 
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same as in the odd order case, but now there is no need of dividing farther 
the arcs 71 and 72, since throughout each region S the sign of the real part 
of pw^ is unaltered. The region Si in which R(pw^ ) si will be considered 
first, that is the integral 



v > 2« J Y X 



np n ' l G{x, s; p n )f(s)dsdp. 



In making use of the expression (3) for G (x , s; p n ) we now alter the last 
column by adding to it each of the columns 1 to /1 multiplied by \yi (s) 
with corresponding index, and each of the columns p. + 1 to n multiplied 
by — \yi (s) with corresponding index. The asymptotic solutions are then 
substituted, and the factor p hi of D is divided into the (1 + i)th row, and 
the factor e^ 6- "' into the jth column. The factor e pw ^ l(b ~ a) is also divided 
out of the (m + l)th column and out of D , thus leaving in the denominator 

A = e pw "( 6_ °> A , 

since «v+i = — «v . Thus G (x , s; p n ) is found to be equal to l/(np™ -1 A) 
multiplied into a determinant of order n + 1 in which the element in the 
first row and jth column is e^*-") [ 1 ] f or j S p. and e l "°' ( - x - b) [ 1 ] for p + 1 
= j = n. The last element is F , given in (29). In the (1 + i)th row and 
jth column the element is [an Wj*] for j < p, [aui w'}'] for p + 1 < j Si n, 
while for j = p, p + 1, and n + 1 the elements are respectively D it Hi, 
and Fi, where D, is as given in (21a), while Hi is the H t given there, multi- 
plied by e l " e ^ b ~ a) , and F and F { are as given in (29). 

The variable s appears only in the last column of this determinant, whence 
f(s) can be multiplied into this last column and the s integration performed. 
Since all terms, save possibly the D's, H's, and F's, are bounded, the expansion 
of the determinant takes the form 



2«J Vi J 
-- f T 



Ff(s)dsdp 

D t Hi fFif(s)ds 

D, Hj fFjf(s)ds 

D r H r fF r f(s)ds 



(a) 



\ i=i 



Ei eP Mi( * _a) 



+ 



—.ft 

2m J i t j=i 



Di 



D„. 



+ £ Ex 

l=n +1 


e W*-»Wp 


,<*-» [ 1 ] 




S t jFif(s)ds 


En , 
-^dp. 


Hj fFif(s)ds 





(b) 



(c) 
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By comparing (a) with the similar integral in the odd order case it is obvious 
that it approaches the limit f(x)/2 as the radius of 71 becomes infinite. 
In (6) the determinants can be written in the form Ep a , where a < 3n. 
Since in the parenthesis the exponents have negative real parts on 71 , (6) falls 
under the type (31), and so approaches zero as the radius of 71 becomes 
infinite. Finally (c) is expanded in terms of the top row of the determinant, 
and the result is treated in a manner similar to that used in the odd order 
case. The expanded form is 



(36) £ie>V->£ = 

i,3 A 



H t JFifds 
H s fFjfds 



i,i A 



D< fFifds 
Dj fFifds 



The formula for J Ft f(s)ds is given in (35) , but a more explicit form for the E 
terms in that expression is now needed. Prom (28) can be obtained directly 
a formula for J Ft f(s)ds. It is 

-—5 f'F i f(s)ds= - £ Wu(Vi) fy 1 (s)f(s)ds 



(37) 



h-X 
1=2 



Z Wu(yj)f 



V}{s)f{s)ds 



y=M+i 



tw li (y j )fyj(s)f(s)d S 

j=l J a 



+1tw ki (y 3 )fy j (s)f(s)<h. 

J=l Ja 



By integrating by parts v times, 



£ 



--pw } {r-a) -— pw/o— a) W 

e-^-Vis)ds --_.[/( T )] +_ s - [/(,)] +7Fl) 



where 



fir) 
I/(r)]-/(r)+^ + 



+ 



f^Hr) 



and [/(<?■)] has the same form. And now this can be applied to each term 
of y } - ( s ) in (26) with the result 



ftp' 



-1 fr .-p»,<T-«) p-pwjitr-a) Tji 

- y i {s)S{s)ds=- zr ^ r [h{r)}+— ^^—[h (<r)] + -^, 

f J* 



Wj J„ - pWj ' • - - ■ pWj 

where the symbol [ h ( r ) ] has the more explicit form 

r it m it \ _i_ ^ l(r) j_ hill. , , ^-i( T ) 

and [ M <r ) ] is the same with the value a substituted in place of r . The 
function h is the function /, but a different notation is introduced because 
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the functions h r that enter are not precisely the derivatives of/. The essential 
thing in the expression is the fact that the functions hr are independent of j . 
When this result is used in (37) and the asymptotic solutions are substituted, 
the result must be of the form (35). Then by comparing the two forms the 
more explicit form for the E terms in (35) is found : 



(38) 



f F t f ( s ) ds m Z ( £ wf t an ] [ h ( a ) ] p*"'-*'" 1 e^"^ 

J 1=2 \j=l 



+ £ ifr[a H ][h(b)]pf 



*!i-*i-l gpW/«|-0) j 



+ ■ 



E 



in which [ an ] is identical with the [ an ] appearing in Z>* (i. e., when the bracket 
symbol is written out in full for any value of m , all corresponding terms as 
well as the E terms are identical). In the same way [an] is identical with 
the [ an ] appearing in Hi . This evaluation of the integrals is now substituted 
in (36), and each determinant there can be written as the sum of several, 
according to the terms of (38). Since from (25) the expressions D t /A and 
( Hi/ A ) are bounded on 71 , the determinants which involve the last column 
E/p are themselves of the form E/p, and that part of (36) is of the form (32), 
and so" converges to zero as the radius of 71 becomes infinite. On the other 
hand the terms involving exponentials other than w„ or w„ +1 reduce to the 
form (31), with one of the constants a* (I 4= 1, h) in place of x. Hence 
these also approach zero. All that is left is then 



E 1 e" w " (x - a) J^^ 
ij A 



Hi Ni 
Hi ^ 



i,j A 



D, Ni 
Di Ni 



where 

Ni = E w*« [ a H ] [ h ( a ) ] p*"-^ 1 e" w ^ tt '- a) 

1=2 

1-1 

+ Z <'Vi lau][h(b)] p*"-*- 1 e p<v+i<»<-» . 



By inequality (25)_the expression (2V.-/A) is of the form E/p. So when the 
values of Di and Hi are substituted from (21a) the first two terms of each 
give expressions that reduce to the form (32) and so approach zero. What 
still remains is 



(39) 
where 



£1* 



,pio Jx-a) 



i,J 



An 

A bl 



B a 



+ Bbi 
+ Ba 



+ 



E 2 e*"V +!<*-*> 



B a 
B a 



+ B w 
+ Bbi 



Aai s £ [ai,]« p "' w (p«;) Mi 



UK, 
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A M m £ [«„•] e'V 4 — > ( - pwrf* h < ( - w„)\ 

!=2 

B« e £ [ 5 K ] [ A ( 6 ) ] p*"-^- 1 e > V s -°<> ( - Wjl ) a « . 



These can be expanded into the sum of a large number of determinants. 
Such determinants as have the same exponential factor in both columns 
are zero, for example 

[ Ob ] e> V 6 -*") ( - p W „ ) A «-^ ( - w „ )*• w h /[a u ][h(b)] p^" 1 e > V 6 -«*> 

[ 5y ] e^*"") ( - pw„ )"'-*> ( - w M )*> w h « [ an ] [ h ( 6 ) ] p V-*M e P>(»-« 2 ) 

for when the common factors 

^<-o|* (t)1( -„„,,.(„,)«.p«.-P«i 

(it is to be noted that [ h ( b ) ] is independent of the subscript) are removed, 

we have 

[an] [au]\ 

which is obviously zero, since the two terms in the same row of the deter- 
minant are identical. 

Since a 2 — a and b — a,k-i are the shortest of the intervals ending at a or b , 
the terms that approach zero most slowly are e<>V ai| - a > and e pw ^ t ~ ak - l) , and 
each of the other determinants in (39) can be written 

Ep 2v ~ 1 e'"V aii '~ a > gP^M^ - "*- 1 ) . 

This divided by A can, by inequality (25), be written E/p, whence (39) is 
of the form (32), and so approaches zero as the radius of 71 becomes infinite. 
To sum up, 

I* = */(*) + e 

for the radius of 71 sufficiently large. 

The treatment of the integral in S 2 , that is over y 2 is, step for step, parallel 
to the above, and the result is obtained that for the radius of y 2 sufficiently 
large 

Iy, = if (*) + «• 

By combining the two results, the convergence of (1) is established in the even 
order case. 



